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Abstract 

We consider the azimuthal correlation of the final-state particles in charged 
weak-current processes. This correlation provides a test of perturbative quan- 
tum chromodynamics (QCD). The azimuthal asymmetry is large in the semi- 
inclusive processes in which we identify a final-state hadron, say, a charged 
pion compared to that in the inclusive processes in which we do not iden- 
tify final-state particles and use only the calorimetric information. In semi- 
inclusive processes the azimuthal asymmetry is more conspicuous when the 
incident lepton is an antineutrino or a positron than when the incident lepton 
is a neutrino or an electron. We analyze all the possible charged weak-current 
processes and study the quantitative aspects of each process. We also compare 
this result to the ep scattering with a photon exchange. 
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I. INTRODUCTION 



The Q CD-improved parton model has shown a great success in describing high-energy 
processes such as deep-inelastic leptoproduction. In the parton model we can express the 
cross section as a convolution of three factors: the parton-lepton hard-scattering cross sec- 
tion, the distribution function describing the partons in the initial state and the fragmenta- 
tion functions describing the distribution of final-state hadrons from the scattered parton. 
The hard-scattering cross section at parton level can be calculated at any given order in 
perturbative QCD. The distribution functions and fragmentation functions themselves can- 
not be calculated perturbatively but the evolution of these functions can be calculated using 
perturbation theory. 

The azimuthal correlations provide a clean test of perturbative QCD since these correla- 
tions occur at higher orders in perturbative QCD. Georgi and Politzer |J proposed the az- 
imuthal angular dependence of the hadrons in the semi-inclusive processes £+p — > £'+h+X, 
where £, £' are leptons, h is a detected hadron. Cahn included the contribution to the az- 
imuthal angular dependence from the intrinsic transverse momentum of the partons bound 
inside the proton. Berger || considered the final-state interaction producing a pion and 
found that the azimuthal asymmetry due to this final-state interaction is opposite in sign 
to that due to the effects studied by Cahn. The azimuthal asymmetries discussed by Cahn 
and Berger are due to nonperturbative effects. These effects were analyzed at low transverse 
momentum. ||.|5[ 

In the kinematic regime attainable at the ep collider at HERA or in the CCFR experi- 
ments, we expect that perturbative QCD effects will dominate nonperturbative effects. This 
is the motivation for considering the azimuthal correlation of final hadrons in ep scattering at 
HERA and in vp {yp) scattering in CCFR experiments. We consider all the possible charged 
weak-current processes in the perturbative regime. Mendez et al. considered extensively 
the azimuthal correlation in leptoproduction. In our paper we analyze the same processes 
but in different viewpoints and analyses. Especially we direct our focus on the experimental 
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aspects since we can now verify the theoretical results in experiments at HERA or CCFR. 

Chay et al. [0] considered the azimuthal asymmetry in ep scattering with a photon 
exchange. Here we apply a similar analysis used in Ref. to charged weak-current processes 
in lepton-hadron scattering. The result is striking in the sense that the final-state particles 
have a strong azimuthal correlation to the incoming lepton. We will systematically analyze 
the azimuthal asymmetry in this paper. In Sec. II we briefly review the kinematics used 
in lepton-proton scattering. In Sec. Ill we define the quantity (cos0) as a measure of the 
azimuthal correlation and calculate it to order a s using perturbative QCD. In Sec. IV we 
analyze numerically the azimuthal correlation in various processes in which the incoming 
lepton is an electron, a neutrino, a positron or an antineutrino. We also compare the results 
from the semi-inclusive processes in which we identify a final-state hadron, say, a charged 
pion with the results from the inclusive processes in which we use only the calorimetric 
information, that is, the energy and the momentum of each particle (or each jet). In Sec. V 
we discuss the behavior of the azimuthal correlation in each process and the conclusion is 
given in Sec. VI. 

II. CROSS SECTIONS 

Here we briefly review the kinematics in lepton-hadron scattering with charged weak 
currents. Let k\ (k 2 ) be the initial (final) momentum of the incoming (outgoing) lepton, 
Pi (P2) be the target (observed final-state hadron) momentum and p\ (P2) be the incident 
(scattered) parton momentum. At high energy, the hadrons will be produced with momenta 
almost parallel to the virtual VT-boson direction, = k^ — k%. We focus on interactions 
that produce nonzero transverse momentum P2T, perpendicular to the spatial component 
of g M , which we will denote by q. We choose the direction of q to be the negative z axis. We 
can write the differential scattering cross section in terms of the following hadronic variables 

Q 2 = -q\ P T = P 2r , 0, 
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Q 2 Pi-q Pi-P2 m 

and the partonic variables 

_ Xh_ Q 2 _ Zh_ _ pi -p 2 / x 

4 2pi • g f pi-q 

The azimuthal angle of the outgoing hadron is measured with respect to kxr, whose 
direction is chosen to be the positive x axis. If we employ jets instead of hadrons, is the 
azimuthal angle of the jet defined by an appropriate jet algorithm || and all the hadronic 
variables are replaced by the jet variables. 

In the parton model, if we consider the inclusive processes £ + p —>■£' + X, in which £, £' 
are different leptons, the differential cross section is given by 

da 



dxndydzftd 2 PT 



J dxdzd£d 2 pT5(xn — £x)5(zh 



x ^(Pt-pt^CQ 2 



d&i 



dxdydzd 2 pT 

Ix H x I ^ T ^ T ^ x ' ^ ^ dxdydzd 2 pT ^ ^ 

with gPP^ = P^dP^d^. The sum i runs over all types of partons (quarks, antiquarks and 
gluons) inside the proton and d&i is the partonic differential cross section. Fi(x,Q 2 ) is the 
parton distribution function of finding the z-type parton inside the proton with the momen- 
tum fraction x. In Eq. (|3]) we neglect the intrinsic momentum due to the nonperturbative 
effects and we identify the momentum of the final-state hadron (or a jet) with the momen- 
tum of the scattered parton. This approximation is valid if we choose final-state particles 
with large transverse momenta. 

If we consider the semi- inclusive process £ + p — > £' + h + X where h is a detected hadron, 
say, a charged pion, the differential cross section is given by 

da 



dxudydz H d 2 PT 



J dxdzd£dt'd 2 p T 5{xH - £x)6(z H - £'z)5 i2) {P T - f'pr) 
da; 



u dxdydzd 2 px 3 ' ^ ^ 



.. Jx H X Jz H Z J Z 



xF '(t^sto d,( t-«")- (4) 

The sum i, j runs over all types of partons. The partonic cross setion da^ describes the 
partonic semi-inclusive process 

£(ki) + parton i(pi) — > ^'(fe) + parton j(p2) + X. (5) 

Here the exchanged gauge boson is a charged W particle. Fi(x, Q 2 ) is the z-type parton dis- 
tribution function, Dj(z, Q 2 ) is the fragmentation function of the j-type parton to hadronize 
into the observed hadron h with the momentum fraction z. These two types of functions 
depend on factorization scales and for simplicity we put the scale to be Q, a typical scale in 
lepton-hadron scattering. 

In order to obtain hadronic cross sections, we have to calculate partonic cross sections 
using perturbative QCD. At zeroth order in a s , the parton cross section for the scattering 
v + q — > e + q' is given by 

da g _ G 2 F m^\V q/q \ 2 Q 2 1 2 



5(l-x)5(l-z)5 2 (p T ), (6) 



dxdydzd 2 pT vr (Q 2 + m^) 2 y 

where V q / q is the relevant Cabibbo-Kobayashi-Maskawa (CKM) matrix element for the pro- 
cess W + + q — > q'. Gp is the Fermi constant and mw is the mass of the W gauge boson. 
For the scattering of an antiquark with a neutrino, v + q — > e + q', the parton cross section 
is given by 



do- q G 2 F m A w \V qq >\ 2 Q 2 (l-y) 
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-5(l-x)5(l-z)5 2 (p T ). (7) 



dxdydzd 2 pT vr (Q 2 + vn^) 2 y 

The only difference between these cross sections in Eqs. (|6]) and (|7|) is the appearance of the 
factor (1 — y) 2 . This is due to the helicity conservation. In short, when particles with the 
opposite handedness scatter, we have the factor of (1 — y) 2 in front, while it is independent 
of y when particles with the same handedness scatter. The cross sections for other processes 
like e + q (q) — > v + q' (q'), V + q (q) — > e + + q' (cf) and e + + q (q) — > v + q' (g 7 ) can be 



obtained using crossing symmetries. However since the transverse momentum is zero at this 
order, there is no azimuthal correlation at the Born level. 

ki / ki 






(a) 



(b) 






(c) 



(d) 



FIG. 1. Feynman diagrams for charged weak-current processes at order a s . 

To first order in a s , the parton scattering processes develop nonzero pt and nontrivial 
dependence on the azimuthal angle 0. The relevant processes are 



q(pi) + W ± *{q) ^ q'{p 2 ) + g{p 3 ), 



(8) 



q{pi) + W ± *{q) 
q{pi) + W ± *{q) 
g(pi)+W**(q) 



q'ips) + g(p2), 
q'(p2) +g(pa), 
q'(ps) + g(p2), 
?(pa) + q'(ps), 
q(ps) + q'(p2), 



(9) 
(10) 

(11) 
(12) 

(13) 



where g is a gluon, W^* is the virtual W boson and q, q' are quarks. The Feynman diagrams 
for these processes are shown in Fig. [I]. Fig. [I|(a) corresponds to Eq. (|8]) [Eq. flT0|)1 with a 
quark line (an antiquark line) and similarly Fig. |l|(b) corresponds to Eq. (Q) and Eq. (pp. 
Fig. |TJ(c) and Fig. |I](d) correspond to Eq. (|T^) and Eq. flT3] ) respectively. 

Using the Sudakov parametrization we can express P2 in terms of x, y and z as 



P2 



(1 — x)(l — z) + xz pi + zq^ + Pt 



(14) 



where pt = (0, p^, 0) is the transverse momentum with pi ■ p? — q • Pt — 0. For massless 
partons we have 



Similarly we can write 



Pt=IP2t| 2 = -(1-x)(1-^)Q 2 . 

X 



y y 



(15) 



(16) 



with fey = (1 — y)Q /y , where fey is defined in the same way as p?- Therefore we have 



h ■ P2 



and 



h -P2 



2xy 



2xy 



[1 - x)(l - z) + xz(l -y) - k r ■ p T 



(17) 



(1 — x)(l — y){\ — z) + xz — • Pr- 



(18) 



The semi-inclusive parton scattering cross section for charged weak-current is given by 
asGlm^lVqigl 2 yQ 2 



d&ij 



dxdydzd?pj 



2tt 3 



(Q 2 + m 2 w ) 2 ' 
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L, u M%5(p 2 T - -(1 - x)(l - z)Q 2 ), (19) 



where L^ v is the average squared of the leptonic charged current and MjT is the partonic 
tensor for the incoming parton i and the outgoing parton j. V q / q are the CKM matrix 
elements. The products L^M^ for the processes in Eqs. @, (|9|), (|i~0|), (|Tl"D, ([12]) and (p~3|) , 
i.e., ij = qq, qg, qq, qg, gq and gq depend on the types of incoming leptons. For the process 
v + parton i — > e + parton j + X, they are written as 

4 (h ■ Pl ) 2 + (k 2 ■ p 2 ) 2 



L (lV M qq 



j A/ft*" 



3 P1-P3P2- P3 

4 (fci ■ pxf + (fc 2 • P3) 2 
3 P\-piPi- P3 
4(fci -p 2 ) 2 + (fc 2 - Pi) 2 

3 P1-P3P2- P3 

4 (fc 1 -p 3 ) 2 + (A; 2 -p 1 ) 2 

3 P\-P1P2- P3 
l (h-p 3 ) 2 +(k 2 -p 2 ) 2 

2 pi • p 2 Pi • P3 

l (fcl-p 2 ) 2 + (fc 2 -p 3 ) 2 
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(20) 
(21) 
(22) 
(23) 
(24) 
(25) 



Pi • P2 Pi ' P3 

Eqs. ( pOQ and (|2l| ) correspond to the Feynman diagrams with quarks in Fig. |I|(a) and |](b) 
respectively with quarks, Eqs. ([22] ) and (p3|) correspond to the same diagrams with an- 
tiquarks. Eqs. (p4|) and fl25|) correspond to Fig. [l](c) and |TJ(d) respectively. Note that 
Eqs. (pT|), ( |23| ) and (^) are obtained from Eqs. (|20|), ( p2|) and fl2"4|) respectively by switching 
p 2 and P3. And Eq. (^) is obtained from Eq. (^) by switching p\ and p 2 . For the process 
e + parton i — > u + parton j + X, the matrix elements squared are the same as Eqs. fl2"0p— fl2"5D 
except an extra factor of 1/2 taking into account the spin average of the incoming electron. 

With the Eqs. (|20|)-(p5|), we can also obtain L^M^ for other charged weak-current 
processes. For example, for the processes V + parton i — > e + + parton j + X, L^M^ are 
obtained by switching k\ and k 2 in Eqs. (^0|)-([25|). They are written as 

4 (k 2 ■ pi) 2 + (h ■ p 2 f 



Pi • P3 P2 • P3 



= 4 (fa -p,)' + (*,■„)' 
y 3 pi ■ p 2 p 2 • p 3 
L MJ », = 4 (fc 2 -p 2 ) 2 + (fci-pi) 2 
M1/ 99 3 pi • p 3 p 2 • p 3 
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(26) 
(27) 
(28) 



L Mt v = 4 (fc 2 ■ p 3 ) 2 + (fci ■ Pi) 
^ qg 3 pi • p 2 p 2 • P3 



(29) 



« - i (fa-^) a +(fa^) 8 , ( 30) 

yv 2 pi • p 2 pi • p 3 

yy 2 Pl • P2 Pi ' P3 

By the same argument L^M^ for the process e + + parton i — > 77 + parton j + X are the 
same except a factor of 1/2. 



III. AZIMUTHAL ASYMMETRY 

The azimuthal asymmetry can be characterized by the average value of cos0, which 
measures the front-back asymmetry of ~P 2 t along the \s.yr direction. It is defined by 



(cos0) = 





! COS0 


J 




1 



(32) 



where da^ (da^) is the lowest-order (first-order in a s ) hadronic scattering cross section 
defined in Eqs. (|3|) or (Q) and the integration over P T , 0, xh, V and zh is implied. When 
we impose a nonzero transverse momentum cutoff, Eq. ( |32]) receives contributions only from 
da^ both in the numerator and in the denominator. Note that the zeroth-order cross section 
is proportional to 5 (Ft)- Therefore with the nonzero transverse momentum cutoff at order 
a s in perturbation theory, the quantity (cos0) is independent of a s . 

In fact the azimuthal asymmetry can occur at the Born level if we include the intrinsic 
transverse momentum due to the confinement of partons inside a proton and the fragemen- 
tation process for partons into hadrons PJ^,0]. However the size of the intrinsic transverse 
momentum due to nonperturbative effects is of the order of a few hundred MeV. There- 
fore if we make the transverse momentum cutoff p c large enough (> 2 GeV) and choose 
hadrons with the transverse momenta larger than p c , we expect that the contributions from 
the intrinsic transverse momentum from the Born-level processes are negligible compared to 
those from a^ 1 ' . In other words the intrinsic transverse momenta of the partons simply can- 
not produce hadrons with transverse momenta larger than p c and the effects from intrinsic 
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transverse momenta are suppressed. Therefore, for p c larger than 2 GeV, (cos0) is given by, 
to a good approximation, 

COS0 



(COS 0) 



da® 



(33) 



In the following analysis we consider (cos 0) as a function of the transverse momentum cutoff 

Pc- 

We first consider the azimuthal asymmetry in the inclusive process u + p — > e + X, where 
X denotes any hadron. The numerator in Eq. (^) can be written as 

da 



where 



da^ COS0 = cos « 



8a s G 2 F m^r 



dxudydz}id 2 PT 
Q 2 l 



X 



3tT 2 (Q2 + m 2 w) 2 y 
/•l (fx 

/ —(A v + B v + C v + D v + E v + F u ), 

Jxif X 



A v 



(1 — y)xz 
\(l-x)(l-z) 



(1 — y)(l — x)(l — z) + xz 



x (|K,| 2 F d (^,Q 2 ) + |K s | 2 F s (^,Q 2 



x 



X 



(1 - y)x(l - z) 



(1 — y)(l — x)z + x(l — z) 



x 



\ (l-a:)z 

(\V ud \ 2 F d (^,Q 2 ) + \V cs \ 2 F s (^,Q 2 ^ 
x x 



a, 



(1 — y)xz 
\{l-x){l-z) 



(1 — x)(l — z) + (1 — y)xz 



x fi^i 2 ^— ,g 2 ) + iK s | 2 ^(— ,q 2 ; 



x 



x 



(1 - y)x(l - z) 



(1 — x)z + (1 — y)x(l — z) 



V \ (l~x)z 
x (\V ud \ 2 F^,Q 2 ) + \V cs \ 2 F- c (^,Q 2 )), 

v X X J 



E„=-{l-2x) 



(1 -y)x(l -x) 



x 



\ z(l-z) 
(\V ud \ 2 + \V cs \ 2 )F g (^,Q 2 : 



z-{l-y){l-z) 



(34) 
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F u = --(l-2x), 



(1 - y)x(l - x) 



x 



>\ z(l-z) 
(\V ud \ 2 +\V cs \ 2 )F g (^,Q 2 



1 - z - (1 - y)z 



The denominator can be written as 
/ da^ = f d 2 P T 



where 



da 



dxudydzud 2 PT 



ActsGpm^ 



Or 



x 



3tt 2 {Qi + mlvfy 
f -(K + Bl + Cl + Dl + El + Fi: 

J xh X 



a: 



B[ 



x 



+ 4(1 ~ v)xz + (1 ~ y)2(1 ~ ~ 
(|K,| 2 F rf (— ,Q 2 ) + |K S | 2 F S (— ,Q 2 )), 

" 1+ ,f (1 7^ + A{1 - y)< 1 ~ z ) + ( 1 - y)^ 1 - x > 

1 (1 — x)z 

(iK,i 2 F d (— ,g 2 ) + iK s | 2 ^(— ,g 2 )), 

v X X J 

C '» = [t 1 ~ V ? (1- + X )(1- Z ) + 4(1 _ V)XZ + (1 " " Z \ 



X 



(|K d | 2 i^(— ,Q 2 ) + |^| 2 ^(— ,Q 

v X X 



(1 - y) 2 J + 4(1 - y)x{\ - z) + (1 - x)z 

(1 — x)z 



X 



E' 



(\V ud \ 2 F^,Q 2 ) + \V cs \ 2 F- c (^,Q'' 

v x x 

l[^+^-vni- ^t{i-z) 

{\V ud \ 2 + \V cs \ 2 )F g (^,Q 2 ), 



3r 



(i-z) 2 + (i-y) 



2^2 



x 2 + (1 



X) 



(\V ud \ 2 + \V cs \ 2 )F g C-f,Q 2 ). 



z(l- 



z 



+ 8(1 - y)x(l - x) 



+ 8(1 - y)x{\ - x) 



(35) 



(36) 



(37) 



The above six terms in Eqs. fl35| ) and ([H]) are obtained from the matrix elements in Eqs. fl20|)- 
(P5]) respectively. For the inclusive process e + p — > z/ + X, the corresponding quantities 
are the same except that the quark flavors are switched, u «-> d and c <-> s in the parton 
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distributions functions. There should also be a factor 1/2 from the incoming electron spin 
average. However it appears both in the numerator and in the denominator, hence cancels 
out. 

Now consider the inclusive process V+p — > e + + X. The numerator and the denominator 
in defining (cos 0) as in Eqs. and (R3) are given by 



and 



Aj7 



(1 - y)xz 



(1 — x)(l — z) + (1 — y)xz 



x 



\(l-x)(l-z) 
(\V ud \ 2 F u (^,Q 2 ) + \V cs \ 2 F c ^,Q 2 )), 

K X X ' 



#77 



(1 - y)x(l - z) 



(1 — x)z + (1 — y)x{\ — z) 



x 



\ (l~x)z 

X H , 1 t t 12 / X H 



(\V ud \ 2 F u (^,Q 2 ) + \V cs \ 2 F c (^,Q 2 )), 
v x x J 



(1 — y)xz 
\(l-x)(l-z) 



(1 — y)(l — x)(l — z) + xz 



x H 



x l\V ud \ 2 F- d (^,Q>) + \V cs \ 2 F- d (^,Q 2 ) 



x 



X 



Fh 



(1 - y)x(l - z) 



(1 — y)(l — x)z + x(l — z) 



x 



\ (l-x)z 
(\V ud \ 2 F d (^,Q 2 ) + \V cs \ 2 F^,Q 2 ^ 



■~(l-2x) t 



x(l -x){l-y) 



x (iK,r+iK s r)F 9 (^,Q 2 ), 



1 - * - z(l - y) 



(1 - - x) r 



x 



8 \ 2(1 — z) 

(\V ud \ 2 + \V cs \ 2 )F g (^,Q 2 ), 

X 



z-{l-y)(l-z) 



(i-y) s 



1 + arz 



2^2 



^,2 



+ 4(1 - y)xz + (1 - x)(l - z) 



x (\V ud \ 2 F u C— ,Q 2 ) + \V CS \ 2 F C (^,Q 2 * 

X X 

B'v = [(i - y) 2 l+ ,f (1 T* )2 + 4(1 - y)x(i - z) + (1 - *)* 

L (1 — xjz 



x (|K,| 2 F U ( — ,g 2 ) + |K s | 2 F c (^,Q 2 



x 



(38) 
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a 



1 + x 2 z 2 



DL 



(l-aO(l-z) 
l + x 2 (l -zf 



+ 4(1 - y)xz + (1 - y) 2 (l - x)(l - z) 



x (iK,r^(— ,q 2 ) + ik s | 2 ^(— ,g 2 ; 



+ 4(1 - y)z(l - z) + (1 - yf(l - X )z 



X 



(1 — x)z 

(\V ud \ 2 F- d (^,Q 2 ) + \V cs \ 2 F^,Q 2 

((1 - zf + (1 - y)V) g2 + (1 ~f )2 + 8(1 - y)x(l - x) 



z(l-z) 



(\V ud \ 2 + \V c f)F 9 (^,Q^ 
' x 



F^=l[{z 2 +{i~ y ni-zf) 

x (\V ud \ 2 + \V cs \ 2 )F g (^ ,Q 2 ). 



x 2 + (1 — x) 



+ 8(1 - y)x(l - x) 



(39) 



For the process e + + p — > V + X, the corresponding quantities are the same as in Eqs. fl38l) 
and (|39|) except the switch of the quark flavors u «-> d and c <-> s in the parton distribution 
functions. 

We can express (cos <fi) using Eq. (|4|) in the semi-inclusive processes in which we identify 
a final-state charged pion. For the process v + p — > e + 7r + X, the numerator can be written 



as 



da^ cos</> 



Q 2 



3vr 2 (Q 2 + m 2 w ) 2 y 

f 1 dx f 1 dz . , , . . 

/ — / — (a„ + b v + + d v + e„ + /„), 



(40) 



and the denominator can be written as 



da (i) = ^a s G 2 F m^ 



Q 2 



x 



3vr 2 (Q 2 + m 2 w ) 2 y 
f - f — « + K + c'u + d' u + e' v + f v ). 



(41) 



The quantities introduced in Eqs. (|4*0| ) and fl4T|) are given as follows: 



(1 - y)xz 



(1 — y)(l — x)(l — z) + 



x 



^(l-s)(l-z) 
(\V ud \ 2 F d (^, Q 2 )DlC~, Q 2 ) + \V CS \ 2 F S (^, Q 2 )D:(^, Q 

v X Z X z 
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z 



c„ = - 



|2 p /%H ^2 



(1 — y)xz 
\(l-x)(l-z) 



(1 — y)(l — x)z + x(l — z) 



x (iK d | 2 ^(— , g 2 ) + ik s | 2 f s (^, g 2 )K(-, g 2 ), 



(1 -x)(l -2) + (1 -y)xz 



x (|K,| 2 ^(-, g 2 )^(-, g 2 ) + \V CS \ 2 F- C (^, Q 2 )D^ Q 2 )) , 

v X Z X Z ' 



(1 - y)x(l - z) 



(l-x)z + (l-y)x(l- z) 



x 



\ (l-x)z 

(\v ud \ 2 F^, g 2 ) + \v cs \ 2 f,(^, q 2 ))d^, g 2 ), 

v X X ' y z 



e „ = -(l-2x). 



(1 - - x) 



z-(l-y)(l-z) 



x (iK.r^^g^ + i^r^-^ 2 ))^— ,g 2 ), 



2j? 



X 



f v = - -{l-2x) t 



(1 - y)x(l - x) 



\-z-(\-y)z 



\ z(l-z) 

x (iK,i 2 ^f g 2 ) + \v cs \ 2 du-, q 2 ))f 9 C— , g 2 ), 



(42) 



where DJ(zh/z,Q 2 ) is the fragmentation function for the i-type parton to fragment into a 
charged pion. 

The quantities in the denominator are given by 

1 + x 2 z 2 



(1 -*)(!-*) 



+ 4(1 - y)xz + (1 - y) 2 (l - x)(l - z) 



b' 



x (\V ud \ 2 F d (^-, Q 2 )DlC-i Q 2 ) + \V CS \ 2 F S (^, Q 2 )D:C-, Q 2 )) , 

v X Z X z ' 

1 + x 2 (l - zf 



X 



+ 4(1 - y)x(l -z) + (l- y) 2 {\ - x)z 
(\V ud \ 2 F d (^, Q 2 ) + \V CS \ 2 F S (^-, Q 2 ))D^, Q 2 ), 



(l-x)z 

X 

1 + x 2 z 2 



(1-1/) 



(l-x)(l-z) 



+ 4(1 - y)xz + (1 - x)(l - z) 



X 



(\V ud \ 2 F^, Q 2 )D^, Q 2 ) + |V CS | 2 F,(^, q 2 )d^, g 2 )) , 



d' 



x " z 

2/1 ,A2 



2 l + x 2 (l-2; 

(1 - y) a-,), 



+ 4(1 - y)x(l - z) + (l- x)z 

|2i7 1 Xh i^2\\^( Zh /02n 



x (|K d | 2 ^(— , g 2 ) + \v cs \ 2 f- c (^ , q 2 ))d-(^, g 2 ), 

v X X ' y z 
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(z 2 + (1 - y)\l - zf) x2+ ( [ l *>" + 8(1 - y)x(l - x) 
v 7 z(l — z) 



x i IK rf | 2 ^(f , Q 2 ) + IV^D^, Q 2 ))F a (^, g 2 ) 



K = ? ((i-^) 2 + (i-y)V 



(1 



x 



(ivy 2 £f ft Q 2 ) + |K S | 2 ^, q 2 ))f 9 { ^, g 2 ). 



z(l-z) 

z 



+ 8(1 -y)z(l-z) 



(43) 



For the process e + p — > I/ + 7T + X, the corresponding quantities are the same as 
in Eqs. ( f42|) and (|43|) except that the quark flavor dependence in the parton distribution 
functions and the fragmentation functions should be switched in each SU(2) weak doublet. 
We can also express the corresponding quantities in the processes V + p — > e + + tt + X and 
e + +p^77 + 7r + X accordingly as in inclusive processes. 



IV. NUMERICAL ANALYSIS 

Let us consider how (cos <fi) behaves numerically when the QCD effects at next-to-leading 
order are included. Note that, if we choose particles with nonzero transverse momentum, 
(cos0) is independent of ct s to first order in a s . Furthermore, if we choose the momentum 
cutoff p c large enough, say, larger than 2 GeV, the contribution of the intrinsic transverse 
momentum inside a hadron is negligible. In our analysis we will show the numerical results 
for the final-state particles with p c > 2 GeV so that we neglect nonperturbative effects. 

We show how (cos0) behaves as a function of the transverse momentum cutoff p c in 
inclusive processes. The numerical results for the inclusive processes with different incoming 
leptons are listed in Table 1. For comparison we list the result from the ep scattering 
in which a photon is exchanged. The plot for (coscp) is shown in Fig. |2|. The numerical 
values are obtained by integrating over the ranges 0.05 < Xh < 0.3, 0.2 < y < 0.8 and 
0.3 < zh{= z) < 1.0. We also require that Q > 2 GeV in order for perturbative QCD to be 
valid. We use the Martin- Roberts-Stirling (MRS) (set E) parton distribution functions 0. 

In Fig. ^| we see that (cos 0) approaches zero as p c increases irrespective of the incoming 
leptons. If we change kinematic ranges, not only the numerical values but also the sign 
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change. However the fact that the azimuthal asymmetry tends to be washed out for large p c 

persists. Therefore the test of perturbative QCD using the azimuthal correlation in inclusive 

processes is not feasible until we have better detector resolution. However in semi-inclusive 

processes the situation is completely different. 

TABLE I. (cos (f>) as a function of the transverse momentum cutoff p c for inclusive processes. The last 
column is from the ep scattering with a photon exchange. The integrated regions are 0.05 < xh < 0.3, 
0.2 < y < 0.8 and 0.3 < z H {= z) < 1.0 with Q > 2 GeV. 



p c (GeV) v — > e e — > v 



2.0 


-0.0192 


-0.0235 


3.0 


-0.0100 


-0.0157 


4.0 


-0.00465 


-0.00979 


5.0 


-0.00200 


-0.00605 


6.0 


-0.000687 


-0.00364 


7.0 


-0.00178 


-0.00194 


8.0 


6.51xl0~ 5 


-0.000952 


9.0 


8.53xl0~ 5 


-0.000355 


10.0 


2.42xl0~ 5 


-0.000107 



v — > e + 


e + — > v 


e — > 


e( 7 ) 


-0.0284 


-0.0192 


-0 


.0351 


-0.0160 


-0.00584 


-0 


.0224 


-0.00852 


0.000910 


-0 


.0145 


-0.00469 


0.00325 


-0 


.00973 


-0.00247 


0.00358 


-0 


.00632 


-0.00116 


0.00277 


-0 


.00401 


-0.000525 


0.00168 


-0 


.00239 


-0.000200 


0.000784 


-0 


.00135 


-6.98xl0- 5 


0.000218 


-0 


.000673 
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FIG. 2. (cos(/>) versus p c in inclusive processes. The leptons listed are the incoming leptons for charged 
weak-current processes. The last one with 7 is from the ep scattering with a photon exchange. 

In the semi-inclusive processes in which we tag a final-state charged pion, we use analytic 
fragmentation functions for simplicity. This is in contrast with studies using Monte Carlo 
simulation for the hadronization process [[HJ. In our numerical analysis we use Sehgal's 



parametrization ||11|| . Sehgal's parametrization for the quark fragmentation functions to 
pions is given by 



D](z) = -(0.05 + 1.05(1 - zf 



(44) 



for j = u, d, u, d and Dj(z) = for other quarks. The gluon fragmentation function to pions 
is given by 



2 2 

D*{z) = -0.1 - 2.1z + — + 4.2 log z. 



(45) 



Note that the gluon fragmentation function is "softer" than the quark fragmentation func- 
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tions, that is, Dg(z) < DJ(z) for z > 0.21. This functional form for the gluon is obtained 
by assuming that the gluon first breaks up into a quark-antiquark pair, and then the quarks 
fragment into the observed hadrons. At large z, the hadrons mainly come from quark frag- 
mentation. For the sake of simplicity, we also nelgect the QCD-induced scale dependence of 
these fragmentation functions. The variation of the fragmentation function due to the scale 
dependence largely cancels out in the ratio defining (cos0). 

Since Q 2 = 2MEiXhV, where M is the proton mass, Ei is the energy of the incoming 
lepton in the proton rest frame, when we integrate over % and y, the strong coupling 
constant a s (Q 2 ) should also be included in the integrands in the definition of (cos0). The 
running coupling constant a s has the Q dependence as 

12% 

a ^ Q2) = (33-2n / )ln(QVA 2 )' (46) 

where rif is the number of quark flavors whose masses are below Q. However the inclusion of 
a s {Q 2 ) in the integrand is numerically negligible since it appears both in the numerator and 
in the denominator. Therefore in our analysis we do not include a s (Q 2 ) in the integrands. 
The numerical error in neglecting the variation of a s with respect to Q is less than a few 
percent. 

TABLE II. (cos 4>) as a function of the transverse momentum cutoff p c for the semi-inclusive processes 
with a final-state charged pion. The last column is from the ep scattering with a photon exchange. The 
kinematic range is 0.05 < x H < 0.3, 0.2 < y < 0.8 and 0.3 < z H < 1.0 with Q > 2 GeV. 
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Pc (GeV) 


v — > e 


e — > v 


v — > e + 


e + — > i/ 


e -> e(7) 


2.0 


-0.0515 


-0.0591 


-0.115 


-0.0817 


-0.0832 


3.0 


-0.0443 


-0.0529 


-0.128 


-0.0854 


-0.0805 


4.0 


-0.0399 


-0.0482 


-0.146 


-0.0970 


-0.0783 


5.0 


-0.0364 


-0.0439 


-0.166 


-0.111 


-0.0762 


6.0 


-0.0341 


-0.0405 


-0.187 


-0.127 


-0.0740 


7 n 

f .u 


n no 1 1 
— U.Uoll 


n nocc 
— U.Uooo 


n on /i 
— U.ZU4 


n i /i 1 
— U.141 


n nvon 
— U.U (2U 


8.0 


-0.0289 


-0.0332 


-0.219 


-0.156 


-0.0698 


9.0 


-0.0256 


-0.0292 


-0.224 


-0.163 


-0.0660 


10.0 


-0.0224 


-0.0248 


-0.226 


-0.173 


-0.0634 



0.00 



-0.05 



A 

01 

o 

o 

V 



■0.10 
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FIG. 3. (cos(/>) versus p c in semi-inclusive processes. The leptons listed are the incoming leptons for 
charged weak-current processes. The last one with 7 is from the ep scattering with a photon exchange. 



The numerical results for the semi-inclusive processes are given in Table 2 and the plot is 
shown in Fig. ^. The numerical values are obtained by integrating over the same range as in 
the analysis of inclusive processes, 0.05 < xu < 0.3, 0.2 < y < 0.8 and 0.3 < zh < 1-0 with 
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Q > 2 GeV. The azimuthal correlation in semi-inclusive processes shows a rich structure. 
As p c increases, (cos</>) decreases for the incoming antineutrino or the positron. On the 
other hand, for the incoming neutrino or the electron, it increases and approaches zero. The 
result from the ep scattering with a photon exchange is located between these two cases. 
This behavior will be analyzed in detail in the next section and we compare it to the behavior 
in inclusive processes. 

V. DISCUSSION 

The most interesting feature of our analysis is the behavior of (cos 0) as a function of the 
transverse momentum cutoff p c . Let us compare inclusive and semi-inclusive cases shown 
in Figs. and |3] respectively. In inclusive processes (cos 0) approaches zero as p c increases 
irrespective of the incoming leptons. On the other hand, (cos0) in semi-inclusive processes 
is numerically large compared to that in inclusive processes by an order of magnitude and 
it depends on the incoming leptons. However (cos <ft) remains consistently negative in semi- 
inclusive processes. Negative values of (cos0) mean that the final-state particles tend to be 
emitted to the direction of the incoming lepton. 

We can understand why there is such asymmetry at order a s in the context of color 
coherence at parton level as noted in Ref. f7j. When a quark-antiquark pair is produced in 
a color-singlet state, soft gluons tend to be emitted inside the cone defined by the quark- 
antiquark pair. In our case, we have an incoming quark and an outgoing quark. However 
we can regard the incoming quark as an outgoing antiquark and the pair as a color singlet. 
Therefore the configuration in which the outgoing quark is closer to the incoming lepton and 
a gluon is emitted between the incoming quark and the outgoing quark is more probable. It 
is this configuration that gives negative (cos </>) after boosting to the photon-proton center- 
of-mass frame assuming that we are in a kinematic regime where the observed hadron is 
coming from the fragmentation of the quark. 

In the semi-inclusive processes in which we identify a final-state hadron, for example, 
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a charged pion, note that the gluon fragmentation function is much softer than the quark 
fragmentation functions. That is, the gluon fragmentation function D^(zh/z) decreases 
rapidly as Zh/z — > 1 compared to the quark fragmentation functions. This is clearly seen in 
Sehgal's parametrization of the fragmentation functions. Therefore for large zh {zh > 0.3 
in our numerical result) we effectively pick up the pions which are fragments of quarks. 
This is exactly the situation where color coherence can explain the asymmetry. Of course, 
final-state quarks can be produced from the gluon- W fusion. But in this case (cos 0) can be 
either positive or negative, hence there is a partial cancellation for wide ranges of xh and 
zh- 

In inclusive processes, since there appear no fragmentation functions, both quarks and 
gluons contribute to the asymmetry. But their contributions tend to cancel each other since 
the final-state particles are emitted in the opposite direction. Note the opposite signs in 
the pairs of terms (A u , B v ), (C u , D v ) and (E u , F v ) in Eq. fl35l). However the asymmetry 
can arise depending on the kinematic range. For example, the valence quarks contribute 
dominantly for large xh/x because the valence quark distribution functions Fi(xi{/x,Q 2 ) 
are larger than other distribution functions. If we compare Figs. and [3], the cancellation 
in inclusive processes is illustrated clearly. The magnitudes of (cos 4>) in inclusive processes 
(Fig. are smaller by an order of magnitude than those in semi-inclusive processes (Fig. [|). 

Now let us consider the detailed behavior of (cos0) as p c varies. In evaluating (cos0), 
there are different combinations of parton distribution functions (and fragmentation func- 
tions in semi-inclusive processes) for different incoming leptons. However, since these func- 
tions appear both in the denominator and in the numerator, the main difference results 
from the matrix elements squared for each process. As the matrix elements squared for the 
incoming electron and for the incoming neutrino are proportional to each other, we expect 
that the behavior of (cos</>) from an incoming electron and from an incoming neutrino is 
similar though the magnitudes may be different. This is true for the cases with an incoming 
positron and an incoming antineutrino. This expectation is shown in Fig. |3| for semi-inclusive 
processes. It is not clear in Fig. § for inclusive processes since the magnitudes of (cos 0) are 
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numerically too small to draw any conclusion. 

One interesting feature in Fig. [| is that when the incoming particle is an antineutrino or a 
positron, (cos <ft) is more negative compared to the case of the incoming neutrino or electron, 
(cos <fi) decreases as p c increases for incoming antileptons, while it increases and approaches 
zero for incoming leptons. This behavior results from complicated functions depending on 
x, y, z, xh and zh- Therefore it is difficult to explain the behavior in a simple way. However 
we can explain why (cos0) is more negative for incoming antileptons with large p c . 

In semi-inclusive processes, since we select the hadron with transverse momentum Pt 
larger than the transverse momentum cutoff p c , we have the relation 

(1 - ,)(! - z) = (1 - - ») 4 , (47) 

XZ xz 

The second equality in Eq. fl47|) is obtained by the relation Q 2 = 2xnyMEi. For large 
Pc the phase space is confined to the region with small x, z and large xh, V and zh- In 
this region the ratio Zh/z, which appears in the fragmentation functions, is large, hence 
the contribution of the gluon fragmentation is negligible compared to that of the quark 
(antiquark) fragmentation. In other words b u , d v in Eq. ( f£^ ) and b' u , d' v in Eq. (^) are 
negligible compared to other contributions. Similarly large xh/x, which appears in the 
parton distribution functions, is preferred hence the contribution of the distribution functions 
of sea quarks and gluons is small compared to that of the valence quark distribution functions 
since the valence quark distribution functions are dominant for large xujx. As a result c u , e u , 
f u terms in (f42|) and c' u , e' u , f' v terms in Eq. (|^) are negligible. Therefore a v and a' v dominate 
for large p c . It means that the main contribution to (cos0) comes from the scattering of an 
initial valence quark into a final-state quark, fragmenting to the observed pion. 

Note that, since the parton distribution functions and the fragmentation functions appear 
both in the numerator and in the denominator, (cos0) is mainly affected by the partonic 
scattering cross sections, which are functions of parton variables x, y and z. For small x, 
z and large y, only the first term in a' v in the denominator and the first term in a v in the 
numerator are important. The partonic part of the integrand in the denominator behaves 
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as (xz)~~ l and that in the numerator behaves as — (xz)~ l l 2 {\ — y) 3 ^ 2 . Since the integrand in 
the denominator grows faster than that of the numerator for small x, z and large y, (cos</>) 
in semi-inclusive processes approaches zero for the incoming electron or neutrino for large 
p c , but it remains negative. 

In the case of the incoming antineutrino, 077 and a' v terms are dominant for large p c as in 
the case with the incoming neutrino. But the behavior of these terms are different. Though 
we do not present the forms of a F and a! v here, we can see the dependence of a v and 
on the partonic variables x, y and z in Eqs. ( |38| ) and ( |39| ) for inclusive processes since the 
partonic cross sections are the same. For small x, z and large y, only the third term in the 
denominator survives and it behaves as {xz)~ x . On the other hand, the integrand in the 
numerator behaves 9bS (xz ) 1/1,2 (1 — y) 1 / 2 . Therefore the magnitude of (cos0) is larger than 
that for the incoming electron or neutrino by a factor of (1 — y)^ 1 in the integrand in the 
numerator, hence (cos</>) is more negative than the case of an incoming electron or neutrino. 
In addition, because of this factor (1 — y)~ l , the difference of (cos0) between the incoming 
antineutrino and the incoming positron is larger than that for the incoming electron and and 
the incoming neutrino. It is also interesting to note that the azimuthal asymmetry exhibited 
by a photon exchange in the semi-inclusive ep scattering is intermediate between the two 
cases in which there are leptons or antileptons. 

The behavior of (cos0) in inclusive processes can be explained by the same argument. 
In this case we identify the transverse momentum of the final-state hadron (or a jet) as the 
transverse momentum of the scattered parton. It corresponds to setting zh = z. Therefore 
we select the final-state particle with the momentum cutoff p c satisfying 



Therefore as p c gets large, the integrated phase space is confined to a region with small x, 
large xh, y and intermediate z between and 1. Since the variable xh/x in the parton 
distribution functions is large, the contribution from the gluon distribution function is neg- 



ligible. This means that E and F in Eqs. fl35f) and (|38|) and E', F' in Eqs. (|37|) and (|39D 




(48) 



x 



x 
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can be neglected. Therefore remaining A, B, C and D terms and their primed quantities 
contribute to (cos (j>) . 

As we can see in Eq. fl36|), the integrands in the denominator behave as x~ l whether 
the incoming particle is a neutrino or an antineutrino. In the case of the neutrino, the 
integrand in the numerator from A v , B u terms behaves as x~ l l 2 (l — y) 3 ^ 2 , while it behaves 
as x -1 / 2 (l — y) 1 / 2 from C u , D u terms. These terms are smaller than the integrands in the 
denominator. Furthermore there is a partial cancellation between A v and B v because they 
have opposite signs. This is also true for C v and D v . Therefore (cos0) becomes very small. 
The same argument applies to the case of the incoming antineutrino. 

As p c gets large, the azimuthal asymmetry tends to be washed out in inclusive processes. 
This behavior of (cos <fi) is expected considering the momentum conservation. In our case in 
which there are two outgoing particles in the VT-proton frame, the transverse momentum of 
one particle is balanced by another particle emitted in the opposite direction. Therefore if we 
sum over all the contributions from all the emitted particles, there should be no azimuthal 
asymmetry. The small azimuthal asymmetry, as shown in Fig. arises since we do not 
include all the emitted particles with the given choice of Xjj, y and z H . 

VI. CONCLUSION 

We have extensively analyzed the azimuthal correlation of final-state particles in charged 
weak-current processes. It is a clean test of perturbative QCD if we make the transverse 
momentum cutoff p c larger than, say, 2 GeV. It turns out that the azimuthal asymmetry is 
appreciable in semi- inclusive processes compared to inclusive processes since the asymmetry 
mainly comes from the contribution of an final-state quark due to the soft nature of the gluon 
fragmentation function for large Zh- In inclusive processes we sum over all the contributions 
from quarks (antiquarks) and gluons, and the sum approaches zero as we include a wider 
range of variables due to the momentum conservation. 

In addition the azimuthal asymmetry is more conspicuous for semi-inclusive processes 
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with an incoming antineutrino or a positron. Previously there was an attempt to analyze 
the azimuthal asymmetry at HERA in ep scattering for electroproduction via a photon 
exchange. However since e + p scattering has been performed at HERA, we expect that the 
test of the azimuthal asymmetry is more feasible because the magnitude of (cos 0) is bigger 
in semi-inclusive processes with an incoming positron. In CCFR experiments they consider 
only the inclusive cross section for (u^) + H — > ji + X, where H is the target 
hadron. If they are able to identify a final-state hadron, they will also be able to observe 
the azimuthal correlations in various charged weak-current processes. 

The azimuthal asymmetry in lepton-hadron scattering results from a combination of 
main ideas in the QCD-improved parton model. As mentioned above, the parton model 
states that the hadronic cross section can be separated into three parts: the parton distribu- 
tion functions, the fragmentation functions and the partonic hard scattering cross section. 
Each element contributes to the azimuthal asymmetry. If we make a transverse momentum 
cutoff p c large enough in order for perturbative QCD to be valid, the small- x (large-x^/rc) 
region mainly contributes, hence the contribution from valence quarks is dominant. At the 
same time, large p c implies that the small- z (large- z H /z) region mainly contributes to the 
asymmetry. This means that quark or antiquark fragmentation functions contribute dom- 
inantly. The detailed behavior of (cos0) depends on the hard scattering cross section at 
parton level. Therefore the experimental analysis of the azimuthal asymmetry tests the very 
basic ideas in the QCD-improved parton model. 
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